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Abstract 

We calculate quantum corrections to the mass of noncommutative </> 4 
kink in (1 + 1) dimensions for intermediate and large values of the noncom- 
mutativity parameter 6. All one- loop divergences are removed by a mass 
renormalization (which is different from the one required in the topolog- 
ically trivial sector). For large 9 quantum corrections to the mass grow 
linearly with signaling about possible break down of the perturbative 
expansion. 



1 Introduction 

The problem of how to compute quantum corrections to the mass of (1 + 1) di- 
mensional solitons was first addressed by Dashen, Hasslacher and Neveu [1] in 
1974. This work was followed by a decade of activities related to quantum prop- 
erties of solitons, see [2,3]. The subject was re-opened in 1997 in the context of 
super symmetric solitons [4-7] (see also reviews [8,9]). About the same time many 
methods of calculations appeared, the ones based on the heat kernel expansion 
are most important for us [10-14]. 

Previous studies of solitons in noncommutative (NC) theories (cf. [15, 16]) 
avoided the problem of quantum corrections to their mass, except for the work 
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[17], where quantum correction to an NC sine-Gordon soliton were calculated in 
a couple of leading orders of the ^-expansion. This paper also reported certain 
difficulties in calculating the mass shift, which are probably not specific to the 
model considered, but rather common to all time-space NC theories. 

The present paper is devoted to calculation of quantum corrections to the 
mass of NC kink in (1+1) dimensions. We do not use the ^-expansion. On 
the contrary, we concentrate on finite and large values of the noncommutativity 
parameter. As in [17] we define the mass shift (the vacuum energy E) as, roughly 
speaking, half the sum of eigenfrequencies of fluctuations about the soliton. Such 
a definitions is hard to justify rigorously since in time-space NC theories where 
there is no canonical definition of energy. Some arguments in favor of treating 
eigenfrequencies as one particle energies can be found in [18]. Anyway, E is useful 
for the discussion of quantum properties of solitons and renormalization. 

This paper is organized as follows. The next section contain some preliminary 
material on the NC kink and fluctuations about it mostly taken from [19]. In sec. 
|3]we study spectral density of the fluctuation and renormalization of the vacuum 
energy (which is done by adding a mass counterterm to the classical action). 
Section [4] is devoted to numerical calculations of the vacuum energy for several 
values of the NC parameter 9 = 2, 3, 5, 7, 10 (in some natural units). The bound 
state frequencies are evaluated with the help of the WKB method, while the 
contribution of the continuous spectrum is calculated by using an approximating 
square well potential. We find that quantum corrections to the vacuum energy 
grow linearly with 9. 

2 Classical solution and fluctuations 

We take the classical action of NC 4 model in 1 + 1 dimensions in the form 



S = / d x 



where (fl = (f) * <ft. Star denotes the Moyal product 



(1) 



exp ( -&"dft> ) f(x)g(y) 



(2) 



where is a constant skew-symmetric matrix which can be chosen as M1/ = 
29e^ w with e 01 = 1. After splitting the coordinates into time and space, {x^} = 
{t,x}, we have the following useful formulae 

f( x ) * e iujt = e lujt f{x + 6u), e iut * f(x) = e lujt f{x - 9u) (3) 
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Obviously, static solutions of the commutative 4 remain also solutions in the 
NC case. In particular, we have the kink solution 



< i m I mx \ , . 

* w= 7x tanl w (4) 



Classical energy of this soliton is 



2\/2 m 
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Qualitatively the spectrum of fluctuations above the kink solution in NC <j) 4 
theory was studied in [19]. Here we repeat some steps from that paper. Let us 
split = $ + (f with $ being the kink solution (j3J) and ip describing quantum 
fluctuations. The linearized equation of motion reads 

- d 2 ip + d 2 x ip + m 2 ip - A($ *$*<p + $*< / ?*$ + < / ?*$*$) = 0. (6) 

Since the kink (J3J) is static, we can look for the solutions of (jS]) in the form 
ip = e lujt r](x). By substituting this ansatz in §6§ and using the relations ([3]) we 
obtain the following equation for rj 

cu 2 r] + d 2 x ri + m 2 r] - A($^ + + ^ 2 _)r] = , (7) 

where $±(x) = $(a;±) and x± = x±8u. It is convenient to choose the mass units 
such that 

m = V2. (8) 

Then eq. ([7]) yields 

{-d 2 x + M 2 + U{x-uo))t 1 = uo 2 7 1) (9) 

where 

U(x;lu) = 2(tanh 2 (x+) + tanh(x + ) tanh(x_) + tanh 2 (x_) — 3) (10) 

and the constant part M 2 = 4 is selected in such a way that U(x; u) — > 
exponentially fast for x — * ±oo and a fixed to. 



3 Quantum energy of the fluctuations 

The vacuum energy can be formally defined as a sum of eigenfrequencies of quan- 
tum fluctuations E = | ^co. In NC theories with time-space noncommutativity 
(in particular, in all NC theories in (1 + 1) dimensions) it is hard to define a 
canonical Hamiltonian (due to the presence of an infinite number of time deriva- 
tives) and thus to justify this formula for E rigorously. Here we follow [17] and 
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accept this definition of the vacuum energy (see the comments made in sec. [T] 
above). More precisely, after taking h — 1 we split the vacuum energy as 



E = E B + E C , (11) 

where 

Eb = \Y,vb (12) 

is a finite sum over the bound state frequencies. The contribution of the con- 
tinuous spectrum has to be regularized. In the zeta function regularization it 
reads 

/■oo 

E c = duo uj 1 - 28 p(uj) , (13) 

J M 

where a regularization parameter s has been introduced which should be put 
zero at the end of the calculations. The function p(u) is the spectral density in 
the continuous spectrum. By using the well known relation between the spectral 
density and the phase shift we can write 

Ec = ^ dk(k 2 + M 2 fi- s d k 5(k), (14) 



where to = \Jk 2 + M 2 . (For the sake of completeness we rederive this formula for 
noncommutative case in Appendix [A]). Equivalently, 

1 f°° 

E c = — dcoco^dJiuj). (15) 



2tt 
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This quantity is divergent in the limit s — > 0. To get rid of the divergences 
one has to know the the asymptotic behavior of the spectral density for large 
frequencies. 



3.1 Asymptotic behavior of the spectral density 

Since the potential in ([9]) depends on the frequency u we are dealing with a 
non-linear spectral problem. To analyze the spectral density we use a method 
developed initially in [20-22] and then adapted to NC theories in [18]. 
Let us consider an auxiliary eigenvalue problem 

L(a)ip a!L0 = u 2 ip a!L0 , (16) 

where 

L(a) = -dl + M 2 + U(x;a). (17) 

The functions r] u = ipu,<j solve our initial eigenvalue problem ([§]), but the spectral 
densities are different. Let us denote the density of the auxiliary problem by 
p(a,u). Note, that according to Appendix lAl we are actually working with the 
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densities from which (an infinite) contribution from free massive fields in an infi- 
nite space has been subtracted. By using this density, one can calculate spectral 
functions of the operators L(a) (with a fixed a). In particular, the heat kernel 
for L(a) reads 

K(L(a); t) = Tr(e"* L(<T) - e -*(-^+^ 2 )) = J du e'^ p(a, w). (18) 

Here Tr denotes a trace over the space of square integrable functions on the real 
line. Again, a subtraction of the "free" heat kernel is necessary. Note, that 
the integration in (|T8|) must be extended over the whole spectrum, including the 
bound states. 

In turn, the spectral density can be expressed through the heat kernel by 
means of an inverse Laplace transformation. 

For any fixed real a the operator L(a) is a Laplace type operator with a 
smooth potential. Therefore, the following asymptotic expansion^ is valid as 
t -> +0 

oo 

K{L{a);t)~J2 tn ' 1/2a ^)- (19) 

n=l 

Odd-numbered coefficients vanish since there is no boundary. The coefficient ao 
does not contribute because of the subtraction of a free heat kernel in ( jl~8l) . Two 
leading coefficients read 

a 2 (a) = -(4tt)- 1/2 / dxU{x;a), (20) 

1 



a i 



(a) = (4tt)- 1/2 / dx 



U{x-af + M 2 U{x;a) 



(21) 



In our case, by using ffTUj) we obtain 
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a 2 (a) = -={9c coth(20a) + 1). (22) 

'71 



For large a we have 
where 



and the corrections in eq. ( 1231) are exponentially small. For the future use we 
note 

a 2j0 = -A f dx f$ 2 -^V (25) 



a 2 (cr) = aa 2 ,i + a 2) o + e.s.t. (23) 



AO 4 

02,1 = — F=! «2,0 = —7=, (24) 

'7T a/7T 



7T / \ A 



^or recent reviews on the heat kernel expansion the reader can consult [23] in the commu- 
tative case, and [24] on NC spaces. 
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Here we used explicit form of the kink solution (j4j) and restored the m-dependence 
on dimensional grounds. 

Higher heat kernel coefficients a 2p are integrals of local polynomials con- 
structed from the potential U (x; a) and its' derivatives. One can easily prove 
that for large a 

a 2p (a) = aa 2pjl + a 2pfi + 0(l/a). (26) 

Probably, the corrections above are even exponentially small, but we shall not 
rely on this. 

The spectral density p(o~, uj) taken at coinciding arguments is not the density 
p(uj) of our initial spectral problem ([9]). As demonstrated in [20-22] (see also [18] 
for a discussion in the framework of NC theories) one has to construct another 
density g(a, to) which is related to a heat-kernel like object 

K(<r,t) = (l + 7 LJL)K(L(*);t) (27) 



2at da , 
through the equation 

K(a; t) = J duo g(a, uj) e~ tuj2 . (28) 

Then 

p(uj) = g(uj,uj). (29) 

We do not know any differential or pseudo-differential operator L(o~) such that 
K(<j; t) = Tr(e~ tL ( a >). In any case, such L cannot be a Laplacian on the real line 
with a smooth potential. However, both K and g are well defined, which allows 
us to consider other spectral functions. 

For t -> +0 



K(<7;i)~^a 2n (<7)H+'< 

n=0 



2 

) 



fl2n(o") = a 2n (o-) + —a 2n+2 (a). (30) 
Za 



For large a we have 



K(a; t) ~ t 1/2 ^2,1 + t l/2 ^aa 2jl + a 2>0 + ^4,1 J + 



(31) 



We are going to use this expansion to evaluate the large uj behavior of g(a,u), 
and, after setting o = uj - the large uj behavior of the physical density p(uo). 
As noted in [22,25], the problem is that besides from powers of uj the high 
frequency asymptotics of the spectral density contain oscillating terms which are 
not defined by the heat trace asymptotics. Strictly speaking, the heat kernel 
expansion defines the asymptotic behavior of the so-called Riesz means of the 
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spectral density rather than that of the spectral density itself. This is not precisely 
what we need, but we learn an important lesson: the power-law asymptotics are 
defined by the heat kernel expansion. Oscillating terms are less important anyway 
since corresponding contributions to the vacuum energy are better convergent 
than that from pure powers. 

Let us consider the zeta function corresponding to g(a, uj) 

((v)= [ du{u 2 )- v e {a,u>). (32) 



One should be careful with possible contribution from oj — 0. One should either 
treat such state separately, or add a small positive part to the mass. The details 
are not essential for us since we are interested in the behavior at large uj. 

There is a well know relation between residues of the zeta function and the 
heat kernel coefficients 

a 2k = Res u= ^ k T(u)C(u) =r(~-k) Res u= i_ k {(v) . (33) 

On the other hand, if the spectral density has a contribution behaving like c p uj p 
at large uj, the zeta function has a pole term 

Jn i/-(p + l)/2 

where Q is a large number, and the coefficient of 2 appeared because we have to 
take into account degeneracy of the continuous spectrum. 

Oscillatory terms do not contribute to the poles. Indeed, after analytical 
continuation from large positive v the integral dujuj~ u sm{buj) has no poles on 
the real line. 

Next we compare (1341) with (1331) and (!30|) . (l3Tj) to obtain the following power 
law asymptotics of the spectral density 

g(a, uj) ~ - 1 a 21 - —= uj~ 2 (a 2tl (T + a 2fi + a 41 / (2a)) + . . . (35) 

For the physical spectral density we have 

p(uj) = g(u, uj) ~ --^-=uj~ 2 a 2fi + 0{uj~ 3 ) (36) 

Zy/TT 



Note, that all terms with uj -1 cancel against each otheiQ This formula does not 
contain possible oscillating term which cannot be obtained by these methods. 

One should not be afraid of negative spectral densities. We have subtracted 
the spectral density of a free massive field. What remains can have both signs. 



2 The terms with l/u> in the spectral density lead to linear divergences in the vacuum energy. 
In the zeta-function regularization such divergences are removed automatically, but they may 
cause problems in other regularization schemes. 
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3.2 Renormalization 



Let us now consider renormalization of the vacuum energy. Since only the con- 
tribution from continuous spectrum is divergent, we treat this term only 



E, 



c 



P 



2s 



duju p(oj) 



(37) 



M 



where we introduced a parameter p with the dimension of mass in order to keep 
proper dimension of the vacuum energy Ec independently of the regularization 
parameter s. Next we choose some frequency Q and split the integral into two 
parts. The part from M to Q we leave as it is. In the part from Q to infinity we 
add and subtract the asymptotics p as (^) of the spectral density. 



E, 



c 



:2s 



P 



duj uj 



M 
oo 



POD 

1 - 2s p(uj)+ fi 2s duji 
Jn 



uj 1 2s (p(uj) - Pas ((jj)) 



+p 



2s 



du uj 1 ' 



(38) 



The choice of Q is simply a matter of convenience. All ultraviolet divergences are 
contained in the last term. As p as we take the first term in fl36l) . 



u a 2 . 



-uj 



7T 



(39) 



where we took into account ([24jl . We assumed that oscillating terms in the spec- 
tral densities (which are beyond our control) do not contribute to the divergences. 
This assumption cannot be universally true. However, in the case we consider in 
this paper the subtraction of ( |39l) indeed gives a convergent integral (see below). 
Therefore, the assumption we made is correct as well as the renormalization of 
the vacuum energy which will be done in a moment. 
The last term in fl38l) can be easily evaluated, 



Em 



div 



.2* 



P 



duj ui 1 ' 



1 



n 



~7T1= fl 2,0 7T 
2\Jn 2s 



py 



Near s 
Ediv 



E div behaves as 
1 



The pole term 



^2,0 



-i + lnf 

2s V 



+ 0(s) 



2 

7i 



2s 



+ In 



+ 0(s) 



(40) 



(41) 



pole 



(42) 



1 m 

ITS y/2 

(where we restored the m-dependence by using dimensional considerations) can 
be removed by an infinite renormalization of the mass 

A 



5m 



2.STT 



(43) 
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c.f. ()5])o. This counterterm does not depend on 9 but is 2/3 of the corresponding 
counterterm in the commutative case [10,12]. 

After removing the pole one can lift the regularization taking the limit s —>■ 0. 



1 f Q 

+ J duu{p{u)- Paa {u))-^-\n(^ (44) 

We remind that p(uo) = (2T[)^ 1 d u 5(uj). 

The presence of a free parameter /i reflects the possibility of a finite renor- 
malization. To fix p one needs a normalization condition. In the commutative 
case in (1 + 1) dimensions it is usually required that the tadpoles are cancelled 
by counterterms (the "no-tadpole" condition). This condition is formulated on 
a constant topologically trivial background. In noncommutative theories one- 
loop divergences on a constant background coincide with the divergences in the 
commutative case and thus differ from the divergences in the kink sector. (The 
difference is precisely the 2/3 factor discussed above). Therefore, the no-tadpole 
condition is not suitable for us. In the case of commutative kink, there is the large 
mass subtraction scheme [29] which is equivalent to the no-tadpole condition [12]. 
Although this scheme can be applied even in non-renormalizable theories, it is 
not clear how to implement it in the noncommutative case. This situation is not 
hopeless, but it is more natural to address it together with studying the commu- 
tative limit 6 — > 0. Our numerical scheme (see below) does not work well in this 
limit. Therefore, we postpone the discussion of physically motivated normaliza- 
tion conditions until a future publication. Here, for the sake of simplicity, we 
put 

p = M. (45) 
The values of E for other choices of p differ by a shift —(2/tt) \n(p/M). 



4 Calculation of energy 

Let us start from the first term in the formula (jUJ) which contains the summation 
over the bound state frequencies. Here we shall distinguish the two cases: first, 

3 In principle one can do renormalization directly in the action by using (|25p . The heat 
kernel expansion for NC </> 4 constructed in [26] predicts the same multiplier of 2/3 which 
relates counterterms in commutative and non-commutative cases. On should however keep 
in mind that the results of [26] are valid for background fields which decay rapidly at the 
infinity. This is not the case of the kink solution, which tends to different constants at two 
infinities. Sensitivity of the heat kernel expansion to the asymptotic behavior of background 
fields is a generic feature of NC manifolds. The only exception are the expansions for operators 
containing only left or only right Moyal multiplications. In this latter case universal formulae 
exist as long as corresponding integrals are convergent [27,28]. 
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Figure 1: The effective potential U(x,uj) + M 2 for uj = 2, 9 = 10 and its square 
well approximation. 



for 9lu considerably larger than 1, and second, 9lu comparable or smaller than 1. 
When 9u) ^> 1 we have very good approximation by a square well potential (cf. 
Fig. Hj). The square well potential is, V — for \x\ > I and V = —V < for 
\x\ < I. (In our case Vq = 4). The scattering data for this potential read 

V e- 2ikl (e 2iu!l - e~ 2i " 1 ) 
821 = Sl2 = (k + uj) 2 e- 2i " 1 — (k — u) 2 e 2i " 1 ^ 

Aujke~ 2ikl 

311 ~ 322 ~ (k + u) 2 e~ 2i " 1 — (k — u) 2 e 2i " 1 ^ ' 

with uj = \Jk 2 + Vq. The width of an approximating square-well potential must 
be chosen such that it correctly reproduces the leading asymptotics of the phase 
shift, i.e. through the condition 

-2Z(w)Vb = J dxU{x-u) = -4(2 + 20u;Coth(20u;)). 

As 21 (u) — 9u + c, it immediately follows c = 1, cf. Fig. HI (Note, that in [19] 
the value c = was taken to estimate the number of bound states for large 9. 
To the leading order in 9 both approximations coincide, but the one chosen here 
reproduces the scattering data of U(x; uj) with a better accuracy.) 
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It is well known from quantum mechanics that the bound states of the square 
well potential obey the equations 



tan(6>a; 2 ) 



0, tan _1 (^ 2 ) + 



0. 



(48) 



to 



In the Table I. one can see some bound states calculated with the above formulas 

(SHD for e = io. 

Alternatively, for finding the bound states one can use the well-known WKB 
formula of the first order 



where X\ and %2 are the turning points. It is well-known that WKB approach 
works well for the low laying bound states, i.e. the smaller 9u, the better accuracy 
of the WKB approximation. At the same time, the larger 9u, the better accuracy 
of the square well approximation. As a result, one can see in the Table I, with 
an example of 9 = 10 case, that highly excited bound state frequencies obtained 
from the square well approximation agree very well with their WKB values in the 
regime of large 9uo. Thus the difference between the WKB and SW (square well) 
data is less than fractions of one percent. Therefore we expect that relative error 
of our calculations of energy for large 9 should not exceed one percent. Let us 
also remind that the zero mode ujq = is the same as in the commutative case. 
(This is the translation zero mode). To summarize, for all bound state frequencies 
LU n >i we use the WKB approximation, which is known to be accurate for lower 
eigenstates, and which practically coincides with the frequencies obtained from 
the square well approximation near the barrier. 

On the contrary to the discrete spectrum, when considering integrals over 
continuous part of the spectrum (second and third terms in (JSj)), we start from 
9 = 2, u = 2 and can use the square well approximation. In other words Q = M 
in (|4*4"j) . if we are limited by not small values of 9. For small 9 and small uj (of the 
bound states and in the beginning of the continuous part of the spectrum) the 
the potential does not look like square well, but rather like a modified Poschl- 
Teller potential. The WKB method is certainly justified for that well for 

calculation of bound states, yet is computationally difficult for the continuous 
spectrum. That is why we did not consider the limit of small 9. 

We used the following expression for the phase shift 



(cf. Appendix |AJ. The dominant asymptotic behavior of p{uj) as at large uj is 
given by (1591 above. 

Since we put Q = M the second term in (14"4"1) vanishes. The calculation of 
the third term in (144"I) was done by numerical integration with a p function given 




(49) 




(50) 
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by a square well approximation. The numerical integration was performed by 
Mathematica and the values for different 9 are given in Fig. 1. One can see 
there that the roughly linear dependence on 9 takes place. The approximating 
square well potential differs from the exact potential U(x; uj) only slightly near 
the points x = ±(9uj + 1), and the form of this difference practically does not 
depend on 9. Therefore, especially since we deal with a massive field, we may 
hope that the total error will remain small and bounded independently of 9, so 
that our conclusion about the linear growth of the integral over the continuous 
spectrum will remain true even if a better approximation is used. Besides, as 
we shall see, the contribution of the continuous spectrum is about 1/2 of the 
contribution of the bound states, so that any error in the continuous spectrum is 
less important. Adding the values for the integration over continuous spectrum 
(Fig. 1) to the sum over the bound state frequencies (see Table I, except for the 
case 9 = 7 which is not presented explicitly in order not too overload the Table), 
we get 



E 


= 2.76, 


(9 = 2) 


E 


= 2.95, 


(9 = 3) 


E 


= 4.66, 


(9 = 5) 


E 


= 6.49, 


(9 = 7) 


E 


= 9.16, 


(9 = 10) 



Let us note that as the contribution to the energy from the integral over contin- 
uous spectrum is linear with 9 (see Fig. 1) and the contribution to the bound 
state is linear as well, the final values for E as a function of 9 is linear in 9 (see 
Fig. 2). 

In the previous section we put \i = M. Values of the vacuum energy for a 
different choice of \x are obtained from the one given in (I5TI) by a constant shift, 
E — > E — (2/tt) \n(n/M). We fixed the mass units so that m = y/2. In arbitrary 
units E should be multiplied by m/y/2. 
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Figure 2: The values of the integral 



— p as (cj)) as a function of 9. 
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Figure 3: The values of the E as a function of 6. The angle of the slope of the 
line is: arctana = 0.726. 
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9 = 


2(WKB) 


9 = 


3(WKB) 


= 


h(WKB) 


9 = W(WKB) 




= 0.955 




= 0.790 


^1 


= 0.613 


LU\ — 


433 


W 2 


= 1.263 


w 2 


= 1.056 


Dn 


= 0.826 


LU2 — 





UJ-i 


= 1.496 




= 1.260 


^6 


= 0.995 


/ 1 In 


fl 708 




= 1.687 


id a 


= 1.143 


U) A 


= 1.135 


1 , t 1 

U/4 — 


n ri 1 

U.Oll 




= 1.847 




= 1.579 




= 1.258 


LU5 — 


n qoi 




= 1.970 


Lila 


= 1.712 




= 1.368 


LUq = 


n 

U. <JOO 






ujj 


= 1.828 




= 1.469 


1 1 Ir, — 
LUJ 


± .KJOO 








= 1.930 


(iJq 


= 1.564 


&8 = 


1 1 98 












= 1.650 


IMg — 


1 1 Q4 










^ 1U 


= 1.737 


UlQ = 












^11 


= 1.810 


LO\\ = 


1 314. 












= 1.880 




1 371 












= 1.946 


^13 = 


1 49^. 














UJ14 = 
















Wl5 — 


1.0Z0 














U)i6 = 


1.0(4 














U>n = 


1.0Z1 














Uis = 


1.000 














CU19 = 


1. < uy 














<^20 = 


1. (01 














1 ilr.-t 

0J21 — 


1 793 














UJ22 = 


1.833 














^>22SW : 


= 1.828 














^23 = 


1.872 














&23SW [ 


= 1.869 














UJ24 = 


1.908 














■ 


= 1.909 














^25 = 


1.944 














^25SW ■ 


= 1.948 














^26 = 


1.978 














&26SW ■ 


= 1.985 



Table 1: Values of the bound state frequencies computed by formula ( HHj) (only 
for the highest five states for 9 = 10) and WKB frequencies found by (jUJ). 
The comparison shows that the WKB approximation results differ from the SW 
(square well) approximation only by fractions of one percent. 
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5 Conclusions 



In this work we found that the divergences in the zeta-regularized one-loop vac- 
uum energy of NC kink (defined as a half-sum over egenfrequencies) can be re- 
moved by the mass renormalization. This renormalization is, however, different 
from the one required in the topologically trivial sector. Although the effective 
potential which defines the spectrum of excitations above the kink depends on 
frequencies, for intermediate and large values of the NC parameter 9 the finite 
part can be calculated by a combination of the WKB method and an approxima- 
tion by a square well potential. For large noncommutativity the one-loop vacuum 
energy grows linearly with 9, so that sooner or later it should become larger than 
the classical value thus signalling break-down of the perturbative expansion. 

Our results may be improved and extended in a number of ways. First of 
all, we need a method of calculations which would work for small 9 and a phys- 
ically motivated normalization condition to fix the parameter /i. Having these 
ingredients at hand, one can address the question whether quantum corrections 
to the NC kink are smooth in the limit 9 — > and whether they reproduce the 
commutative result [1] in this limit. 



A Vacuum energy in the zeta regularization 

Here we derive eq. (|T^|) by using the approach of Bordag [10] and making necessary 
modifications due to the noncommutativity. In our exposition we also use Ref. 
[12]. First we introduce a cut-off at large distances by imposing the Dirichlet 
conditions on the perturbations r](—L(k)) = rj(L(k)) = 0. In commutative case, 
when the potential has no dependence on to, it is enough to take L{k) = const. 
In our case the potential does depend on uj. We would like to ensure that x± are 
far away from the boundary for all uj. This can be achieved by taking L(k) = 
9uj + Lq where L is a large positive constant. Later we shall consider the limit 
L — > oo. The frequency dependent boundary condition is the main novelty here 
as compared to previous works. We shall demonstrate that it does not change 
the result. 

Without boundaries for each momentum k there are two independent solutions 
7)i, 7)2 of the wave equation with the asymptotic behavior 

r)i ~ e ikx + s 12 e~ ikx , 7)2 ~ s 2 2e ikx for x -> -oo 

Vi ~ s u e ikx , 7)2 ~ s 2 ie ikx + e~ ikx for x -> oo. (52) 

The potential U is symmetric under the reflection x — > —x. Consequently Su = 
S22, «2i = s i2- For large but finite L the spectrum is discrete and is defined by 
the condition 

f(k) = (( Su + s 21 y kL + e- ikL )(( Sll - s 21 y kL - e~ lkL ) = 0, (53) 



15 



where the bracket with the plus (resp., minus) sign corresponds to a symmetric 
(resp., antisymmetric) solution. 

It is known that if we have discrete spectrum only, the zeta-regularized vacuum 
energy is a sum over the eigenfrequencies, E = \ J2 n w i~ 2s = \ J2 n (^n + M 2 )?~ s . 
The function dk In f(k) has poles at k = k n with unit residues. Therefore, we can 
rewrite regularized Eq in the form of a contour integral, 



^AimV+^m*™- (54) 



The integration contour runs anticlockwise around the real positive semiaxis and 
consists of one branch at k = Re k + ie, a second branch at k = Re k — ie, and a 
small segment — e < Im k < e along the imaginary axis. Along the upper part of 
the contour we keep in f(k) only the terms with exp(—ikL(k)) since exp(ikL(k)) 
vanishes as L — > oo. Along the lower part of the contour we retain exp(ikL(k)). 
The contribution from the third part can be dropped, as in [10,12]. We have 

4 L0) = \[^ + M^-^(iikL(k) + ln(4 - s\ 2 )). (55) 



Next we take into account 



^-4 = e M(fc) , (56) 



where 5{k) is the phase shift, and subtract the contribution from free fields of 
mass M satisfying the same boundary conditions (i.e., the expression ( 1551) with 
S(k) = 0). After taking the limit L — > oo we obtain eq. (Tl4|) for the vacuum 
energy. 

Note, that taking the boundaries into account explicitly is essential in su- 
persymmetric theories (where the boundary conditions must be supersymmet- 
ric) [12]. 
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